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Abstract — The concept of plug-in electric vehicles (PEV) are 
gaining increasing popularity in recent years, due to the growing 
societal awareness of reducing greenhouse gas (GHG) emissions, 
and gaining independence on foreign oil or petroleum. Large- 
scale deployment of PEVs currently faces many challenges. One 
particular concern is that the PEV charging can potentially cause 
significant impacts on the existing power distribution system, due 
to the increase in peak load. As such, this work tries to mitigate 
the impacts of PEV charging by proposing a decentralized smart 
PEV charging algorithm to minimize the distribution system load 
variance, so that a 'flat' total load profile can be obtained. The 
charging algorithm is myopic, in that it controls the PEV charging 
processes in each time slot based entirely on the current power 
system states, without knowledge about future system dynamics. 
We provide theoretical guarantees on the asymptotic optimality 
of the proposed charging algorithm. Thus, compared to other 
forecast based smart charging approaches in the literature, the 
charging algorithm not only achieves optimality asymptotically in 
an on-line, and decentralized manner, but also is robust against 
various uncertainties in the power system, such as random PEV 
driving patterns and distributed generation (DG) with highly 
intermittent renewable energy sources. 

Index Terms — Distribution systems, smart charging, on-line 
algorithm, plug-in electric vehicle, minimum load variance, smart 
grids. 

Nomenclature 

Ai (n) Energy consumption of PEV i at time slot n. 

yl™ Maximum energy consumption of PEV i dur- 

ing a time slot. 

Ci Charging reference offset for PEV i. 

f(d) Daily average charging cost for day d. 

jmax Maximum charging cost in a time slot. 

N Number of customers in the system. 

Pi (n) Charging power of PEV i at time slot n. 

pmax Maximum charging power of PEV i. 

<5 n et(^) Total net base load observed by the substation 
at time slot n, which may include the output 
of distributed generators as negative load. 

Snet* Maximum total net base load in a time slot. 

At Length of a time slot. 

T Number of time slots for one day. 

Ui(n) Energy queue length of PEV i at time slot n. 

U ref (n) PEV charging reference signal set by the ag- 
gregator at time slot n. 

ft Weight of the charging cost. 

rji Charging efficiency of PEV i. 
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/j,(d) Daily average of the total net load for day d. 

Xi(n) Indicator function of whether PEV i is avail- 

able for charging at time slot n. 

I. Introduction 

THE growing societal awareness of environmental issues, 
as well as ongoing concerns about reducing the depen- 
dence on foreign oil or petroleum, have made the concept 
of plug-in electric vehicles (PEV) very popular during the 
past few years [1]. PEVs, such as electric vehicles (EV) and 
plug-in hybrid electric vehicles (PHEV), can contribute to 
resolving these energy security issues by reducing the green- 
house gas (GHG) emissions and petroleum consumption in the 
transportation sector. Currently, large-scale implementation of 
PEVs in the near future is being planned. For example, the 
present US administration has planned 1 million PHEVs by 
2015. Realizing the critical role of PEV in the auto industry, 
many automakers, such as Toyota, Nissan, GM, BYD, Fisker 
and Tesla are planning to produce dozens of types of PEVs 
(TJ, with many coming out starting by 2012. 

It has been widely recognized that high penetration level 
of PEVs will cause significant impacts on the existing power 
system, in particular at the distribution level 0, 0, 0, 
13, 0, 0, 0, 0. Without proper coordination, it will 
be very likely that most of these PEVs will start charging 
during the overall peak load period 0, causing severe branch 
congestions and voltage problems. Some studies 0, have 
shown that the existing distribution system infrastructure may 
only support a very low PEV penetration level (such as 10%) 
without grid operation procedure changes or additional grid 
infrastructure investments. On the other hand, studies have also 
shown the promising result of mitigating the impact of PEV 
charging by coordinated charging, which can effectively shifts 
the PEV load to the off-peak period. For example, it has been 
shown that coordinated charging can achieve around 50% PEV 
penetration level in certain existing distribution systems 0. 
Thus, it is crucial to design effective coordinated PEV charging 
algorithms for large-scale deployment of PEVs in the current 
power system, in order to achieve efficient grid operation, as 
well as bypassing or deferring the costly grid infrastructure 
investment. In other words, the existing 'cyber' infrastructure 
may be used to resume the 'physical' infrastructure in the 
power cyber-physical system (CPS). 

This work contributes to the integration of the PEV into 
the power system by proposing an on-line smart charging 
algorithm. The algorithm tries to minimize the distribution 
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system load variance. The reason that such an objective 
function is chosen is as follows. Firstly, the minimum load 
variance objective function can achieve a perfect 'valley- 
filling' charging profile, in the sense that at the optimal solu- 
tion, the total load profile is as flat as it can possibly be 1 10 1, 
ifTD . This implies that the PEV charging load can be efficiently 
'spread' among the off-peak periods, which may help achieve 
a higher PEV integration level in the existing power system, as 
well as lowering the distribution system loss 1121 . as compared 
to the other smart charging algorithms, in particular the ones 
based on electricity price (6), fl3l . In the latter case, it is 
possible that a new 'PEV charging peak' can form during 
the midnight, as many PEVs start charging simultaneously, as 
triggered by low electricity price. Secondly, the minimum load 
variance formulation is convex, and therefore can be solved 
exactly and quickly. Lastly, the proposed objective function is 
very flexible, which can also be used for other applications, 
such as unidirectional vehicle-to-grid (V2G) ancillary services. 
This can be done by 'modulating' the base load according to 
external regulation signals. The details will be discussed later 
in this paper. 

As the main contribution of this work, this paper proposes 
an on-line decentralized algorithm to solve the minimum 
load variance PEV charging problem. The charging algorithm 
is easy to implement. In each time slot, 'binary' charging 
decisions (charge or not charge) are made locally by each 
vehicle after comparing its battery's state of charge (SoC), 
which is the available percentage of the battery capacity, 
to a charging reference signal set by an aggregator, which 
is owned by the utility. The charging reference is carefully 
chosen, based entirely on the current states of the power 
system, such as the SoC values of the plugged-in PEVs, the 
output of distributed generation (DG), the household base 
loads, as well as external regulation signals. In particular, the 
charging reference signal is chosen to greedily optimizes a 
function of the current system state, maximizing a battery 
'energy queue' weighted charging power, penalized by a 
quadratic function of total load in each time slot. Somewhat 
surprisingly, we will show that such a myopic charging algo- 
rithm achieves the same asymptotic performance as compared 
to any optimal scheduling algorithm with perfect one-day 
ahead forecast of all uncertain parameters. In other words, 
our algorithm achieves the same optimality as compared to 
conventional approaches, such as dynamic programming |14|, 
while dramatically reducing the computational complexity. In 
fact, the proposed charging algorithm is closely related to 
the celebrated optimal max-weight policy in the stochastic 
control literature lfT31 . fl6l . Such max-weight type algorithms 
are intimately related to the stochastic sub-gradient algorithm 
1 16 1 and dynamic programming |fT31 , and has found numerous 
successes in diverse areas such as computer networks fTTIl , 
wireless networks [18], [19] and power transmission system 
1 15 1 . In this paper, we extend the application of the max- 
weight algorithm to the PEV charging problem in the power 
distribution system, and prove the optimality results. 

The smart PEV charging is an emerging area of research, 
which has been subject to continuing investigations. While 
there are many algorithms addressing smart PEV charging, 



most of them are solved in a centralized manner Q], ||6l , 
lfl2l . ||20l , II2TI . assuming sufficient accuracy on day-ahead 
predictions about the stochastic dynamics in the power system. 
The real-time, decentralized charging issues are addressed by 
very few works, such as [22], [23], where it is very challenging 
to obtain performance guarantees. Note that the forecast based 
algorithms are vulnerable to the prediction errors. In particular, 
they may face significant performance degradation as the 
current distribution system gradually evolves into the future 
smart distribution system [24) , where large-scale integration 
of distributed generation with intermittent renewable sources 
can result in highly uncertain stochastic dynamics, which are 
hard to predict from one day ahead. On the other hand, the 
proposed PEV charging algorithm in this paper does not suffer 
from such performance loss, as the PEV charging decisions 
are determined completely in real time, where power system 
states can be observed with very good accuracy. Thus, the 
proposed PEV charging algorithm can achieve both optimality 
and robustness in an on-line, decentralized manner. 

The rest of this paper is organized as follows. In Section 
[TT1 we introduce the model and formulate the optimal PEV 
charing problem. Section[In]proposes and analyzes the optimal 
PEV charging algorithm, and Section [IV] demonstrates the 
simulation results. Finally, Section [VI concludes this paper. 

II. Assumptions and Modeling 

In this section, we introduce the system model and formu- 
late the minimum load variance PEV charging problem. 

A. Battery State Model 

A discrete-time system is considered in this paper, where 
the length of each time slot matches the typical sampling and 
operation time scale of the PEV aggregation unit. For example, 
the sampling rate of the distribution system load can be on 
a 15-min time basis [4|. The PEV aggregator is assumed to 
be owned by a utility company, which is concerned about 
potential grid problems with PEV charging. Thus, the goal 
of the PEV aggregator is to achieve on-line, decentralized 
coordination of PEV charging, so as to minimize its impact 
on the distribution system. 

Denote T as the total number of time slots for each day. We 
assume that there are N customers in the power system, and 
denote S net (n) as the total net base load as observed by the 
aggregator during time slot n, which may include the output 
power from distributed generation as negative loads: 

S net {n) = Sbase(n) - S DG (n) (1) 

where Sbase and Sdg are the total base load and total distributed 
generation, respectively. Denote Ui(n) as the 'energy queue 
length' at PEV i at the beginning of time slot n, which is the 
amount of energy that needs to be charged to refill the battery 
of PEV i. That is, 

Ui(n) = (1 - SoC % {n)) x Cap, (2) 

where SoCi and Cap i are the state of charge and capacity of 
PEV i's battery, respectively. Thus, the queueing dynamics of 
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the Ui(n) can be expressed as the following: 

Ui{n) = Ui{n - 1) - rnPi(n)At + A^n) (3) 

In above, -Pj(n) is the charging power of PEV i during time 
slot n, At is the slot length, and rji is the battery charging 
efficiency of the PEV i, which depends on the PEV charger 
and battery types. It is assumed that 



< Pi(n) < i^ max ,l < i < N 



(4) 



where p* nax depends on the charging circuit rating for PEV 
i. For example, for the standard 120-V/16-A wall outlet, the 
maximum charging power is 1.92kW. One technical issue 
is that the 'energy queue length' cannot be negative, which 
implies that 



rnPi(n)At < Ui(n - 1), 1 < i < N 



(5) 



This simply specifies that a battery will not further charge 
once it is full. Ai(n) is the random energy consumption of 
the PEV during time slot n. Bidirectional V2G applications are 
not considered in this paper, and will be addressed in future 
research. In this paper, the energy consumption Ai (n) can be 
only caused by the random PEV driving activity J9) at time 
slot n, during which the vehicle is not able to draw power from 
the distribution system. Such charging availability is formally 
described by Xi( n )' which is an indicator function of whether 
PEV i is plugged into the power grid at time slot n, i.e., 
Xi( n ) = 1 if PEV i is plugged into the power grid at time 
slot n, otherwise Xi(n) = 0. Thus, we have Pi(n) = if 
Xi( n ) = 0' since the PEV can not draw any power from the 
grid when it is not plugged into the power grid. 

B. Minimum Load Variance PEV Charging 

This paper focuses on coordinated PEV charging to mini- 
mize the total load variance as seen by the PEV aggregator. 
The optimal coordinated PEV charging problem can be for- 
mulated as follows: 



^ T N 

s.t. 

P<(n) = 0if Xi (n) = 0,Vi,n 
Ui(T) = Ui(0),l<i<N 



(6) 

(7) 
(8) 
(9) 



In above, the cost function / in © corresponds to the total 
load variance as seen by the PEV aggregator, where 
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T N 

71=1 1=1 



Aj(n) 
rji At 



(10) 



is the average total load during one day. The constraint in 
(O essentially requires that all PEV energy consumptions 
during the day should be met by the PEV charging schedules 
{Pi(n)}. Note that our formulation has the same objective 
function as the minimum load variance formulation in lfT2l . 
with more detailed modeling of the battery states and PEV 
driving processes as in (fj) and ©. Further, note that both 
sequences of the net base load {S net (n)} an d the PEV driving 



variables {xi( n )}> {Ai(n)} are stochastic processes. In the 
literature, such randomness is treated either as deterministic 
1 13 1, or by stochastic programming methods [4], which require 
the knowledge of the joint probability distribution of correlated 
random variables. 

One obstacle in achieving completely myopic, real-time 
PEV charging is that the minimum load variance in ©-Q 
requires knowledge of the average load /i, which is a function 
of PEV loads during the whole day. On the other hand, such 
average load information is not needed to solve ©I-©, as we 
show in the following theorem: 

Theorem 1: Consider the following optimization problem: 



T N 

min /= A £(S net W + P(n)f 

{P,(n)} 1 • 



S.t. 



71 = 1 
0,®,© 



i=l 



(ID 



Any optimal solution of ( fTTT > is also optimal for (©-(O. 

Proof: See Appendix lAl ■ 
Thus, in order to solve the original minimum variance PEV 
charging problem, it is sufficient to solve ( fTTT i instead, which 
does not include the average load /i, and therefore allows easy 
implementation of myopic algorithms. Throughout of the rest 
of the paper, we will directly focus on solving ( TlTT i, bearing in 
mind that the original problem ©-(O is also solved, according 
to Theorem [TJ 

In the literature, most smart charging algorithms (6), 
[10 1, 1 12 1 solve the optimization problem ( fTTT > directly, as- 
suming that all the random variables {S net (n)}, {Ai(n)} and 
{Xi(n)} (or their joint probability distribution) are known. 
Thus, the performance of these algorithms depends crucially 
on the prediction accuracy of these parameters, requiring 
knowledge of automotive driving patterns and base loads. 
While the former is clearly unrealistic, even loads may be 
unpredictable, if the system has a large number of small- 
scale distributed generators, which can use highly intermittent 
renewable energy sources, such as wind and solar. Prediction 
errors can cause substantial performance loss in these charging 
algorithms. As an alternative approach, in this paper, we 
propose a myopic PEV charging algorithm to solve (fTTl i. which 
can achieve the optimal charging cost asymptotically, without 
suffering from day-ahead prediction errors. 

III. Decentralized Smart PEV Charging 

In this section, we describe and analyze the performance of 
the on-line decentralized PEV charging algorithm. 

A. Decentralized PEV Charging 

The algorithm is shown in Algorithm Q] According to the 
algorithm, the charging decisions are made by each vehicle 
locally by comparing its battery 'energy queue length' (or 
equivalently, its SoC) to a charging reference set by the 
aggregator. Fig. Q] illustrates such local charging process at 
a PEV in the power system. In this way, the PEV charging 
processes are properly coordinated, so that the PEVs with 
higher 'charging pressure' Ui(n — 1) (equivalently, lower 
SoC) can be served first. Further, compared to the centralized 
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Algorithm 1 On-line Decentralized PEV Charging 
1: Initialization: At the beginning of each time slot, the 
controller of each plugged-in PEV initializes its charging 
power Pi <— 0, and the aggregator initializes the global 
charging reference as follows: 



' ref 



: )/2 



(12) 



where U™" 1 and L/^ ax are lower and upper bounds, 
respectively, which are selected based on historical data. 
2: Each vehicle controller updates its charging power: 



Pi 



Ui(n - 1) < £W(r?iA*) - Q 
pmax otherwise 



(13) 



where Ci is a properly chosen constant depending on the 
charging service contract of PEV i. 
3: Each vehicle controller submits the intended charging 
power Pi to the aggregator, which sets £/^ ax <— U te f if 



N 



U Kt >2f3{S net (n)+J2 P i. 



(14) 



Otherwise, it sets U£g n 



U Te {. P is a properly chosen 
constant to achieve a trade-off between refilling the bat- 
teries quickly and minimizing the charging cost. 
4: The charging reference C/ re f is broadcast to each vehicle 
controller. If \U£g n - U^ ax \ < e', stop. Otherwise go to 
Step 2. 



charging scheduling approaches, the algorithm is also less 
intrusive, since the PEVs are the decision makers, and are 
only required to submit their charging power {Pi} to the 
aggregator. Note that neither the charging decisions in ( TT3b nor 
the charging reference U K f update is done in an ad hoc fashion. 
In fact, we will show the iterative procedures in Algorithm [TJ 
solve a generalized 'max-weight' optimization problem. 

Theorem 2: In each time slot n, the charging power 
{Pi(n)} computed by Algorithm [TJ solves the following opti- 
mization problem: 



JV N 

m&xJ2(U t {n - 1) + C^rjiPiAt - (3(S mt (n) + J2 P * 

{ Pi ^ i=l i=l 



s.t. < Pi < if"*, 1 < % < N 

Pi^At < U l {n - 1), 1 < i < N 
Pi =0 if Xi[n) = 0, 1 < i < N 



(15) 



For the case that the battery of PEV i is fully charged in 
the middle of a time slot, Pi(n) is interpreted as the 'average 
charging power' Pi(n) = Ui(n — l)/(rjiAt) for that particular 
time slot. 

Proof: See Appendix IB] ■ 
We will continue with the above analysis in Section IIII-CI 
and present performance guarantees of Algorithm Q] Before 
that, we will discuss the intuition behind the algorithm, and 
its key properties. 



Charging Process at Vehicle 1 



Quene Length 
Threshold 

Charging Power (kW) 




Time (h) 



Fig. 1. Illustration of the decentralized charging profile of a typical PEV. 



B. Discussions on the Charging Algorithm 

1) Intuition behind the Algorithm: Note that at the con- 
vergence, the charging reference can be written as U K f = 
2(3(S net (n) + Yli=i Pi)- Thus, intuitively, during the peak 
period, the aggregator will set a high charging reference U le f, 
so that the PEVs with low battery 'energy queue length' 
(correspondingly, high SoC) will not charge. Similarly, for 
the off-peak period, the aggregator will set a low charging 
reference t/ le f, so that the PEVs will 'fill the valley'. Such an 
approach can be easily used to absorb generations with highly 
intermittent renewable sources, by treating them as negative 
loads in S net (n). Note that it is also possible to 'modulate' the 
net base load S 1Kt (n) according to certain external regulation 
signals, in order to achieve V2G ancillary services. Detailed 
discussions on the specific methods, on the other hand, is out 
of the scope of this paper. 

2) Generalized Max-Weight Policy: Anther way to interpret 
the algorithm is to inspect the optimization in (fT5l l. Note 
that the objective function in ( TT3T > has two parts: the 'energy 
queue length weighted charging' term (Ui(n—l) + Ci)rjiPiAt, 
and the 'charging cost' term (S ne( (n) + 2~2tLi Pi) 2 * which 
is weighted by parameter (3. The first term is related to the 
'stability' requirement of the energy queues {Ui(n)}, which 
encourages the PEVs to charge at high power, in particular 
the ones with low SoC. On the other hand, the PEV charging 
is also penalized by the second quadratic term, which is the 
instantaneous charging cost at time slot n, so as to prevent 
grid problems due to aggressive PEV charging. Finally, by 
adjusting the parameter f3, one can achieve a trade-off between 
quickly charging the PEV batteries (maximizing the first term), 
and achieving minimum charging cost (minimizing the second 
term). The effect of f3 will be illustrated in case studies. 
Note that (3 needs to be chosen by the aggregator carefully, 
according to system specifications and historical data. 

3) On-Off Charging: Algorithm [TJ has an interesting 'on- 
off charging property. That is, each PEV either charges at 
the maximum power, or do not charge at all. This can be seen 
from the algorithm specification in d!3t . as well as in Fig. [TJ 
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Fig. 2. The communication structure of the decentralized PEV charging 
algorithm. 



Thus, compared to the smart charging algorithms, which adjust 
charging powers continuously, the on-line decentralized PEV 
charging simplifies PEV charging circuit design. Similar type 
'on-off charging profiles are also considered in 1251 . 

4) Communication Structure: The communication structure 
for the algorithm is shown in Fig. |2] According to this 
structure, the vehicle controllers directly communicate with 
the aggregator, perhaps through extra intermediate levels of 
aggregation (6). There is no direct communication between 
the PEVs. Further, the aggregator only needs a scalar, namely 
the aggregated PEV load (X^i=i ^i)> an d broadcast back the 
charging reference U K f to the PEVs. Therefore, the proposed 
scheme requires low communication capability, and is easy for 
real-time implementation. 

C. Performance Guarantees 

This section continues with the analysis of Theorem [2] 
and demonstrates the asymptotic optimality of the proposed 
PEV charging algorithm, by comparing it against the optimal 
solution of ( fTTT ) over a long time period. As a minor technical 
assumption, we assume that the PEV energy consumptions 
{Ai(n)} are strictly feasible, in the sense that there is a 
charging schedule, which satisfies (|4), (|8), and that 

T T 

Pi(n)rnAt > Vii(n) +sT,l < i < N (16) 



71=1 



for a small constant e > 0. This simply means that the duration 
that a PEV is plugged into the distribution system is 'more 
than just enough' to refill its daily energy consumption. We 
are interested in the performance of the algorithm over a period 
of D days, and denote f*(d) as the optimal cost of (fTTb for 
the d-th day. We have the following theorem: 

Theorem 3: Algorithm Q] achieves the following asymptotic 
average charging cost: 

^ DT N 

limmp — + S net (n)) 2 < 
D ^°° DT ,7=! 



and the following average energy queue lengths: 

Inn sup jff;2^2^ U A n ) < KZ 



D-s-oo 



71=1 7=1 



Bi , (3f n 



N 



d=l 



7=1 

where the constants are defined as follows: J max = 
(^t ax + E7li^ max ) 2 -/i = E^i(Ci + Ki)K it B 2 = 
Y^ n =i eK i< B 3 = J2n=i K i> and the constant K t = 
max(^ nax ,77 lJ Ff lax At). 

Proof: See Appendix ICl ■ 

From the above theorem, it is clear that, by properly 
adjusting the parameter j3, the long term average charging 
cost can be made arbitrarily close to the optimal, with an 
optimality gap on an order of 0(l/f3), as shown in ( fT71 ), 
while the average queue lengths grow as 0((3), as shown in 
(list . The asymptotic bound provided by the above analysis 
is not tight in general. On the other hand, it is possible to 
provide tighter guarantees, by assuming specific probabilistic 
distributions on the stochastic processes, such as PEV driving 
patterns, base load, and renewable generation. Such a model 
dependent analysis is out of the scope of this paper, and will 
be addressed in future research. 

While this section showed the theoretical optimality of the 
decentralized PEV charging algorithm, in the next section, 
we demonstrate its performance in simulation to illustrate this 
point. 

IV. Case Study 

In this section, we test the performance of the proposed 
algorithm, and compare it against other prediction based 
algorithms. The simulations are done on IEEE 37-bus system 
and IEEE 123-bus system |26|, with MATLAB on a Xeon 
X3460 CPU with 8GB of RAM. 



A. Simulation Setup 

1 ) Base Load: The test systems are assumed to be residen- 
tial distribution systems. Both systems use the same household 
base load curve, which is chosen according to the typical 
southern California residential load from the SCE website [27|, 
with proper scaling to match the national average household 
load of 1.3kW [28|. Fig. [3] shows the daily curve of total base 
load used for the simulation with the IEEE 37-bus system. 
The IEEE 123-bus system uses a properly scaled curve with 
the same shape. For day-ahead prediction based PEV charging 
algorithms, we assume that they use forecast profiles such as 
the one in Fig. [3] Note that the shaded region corresponds to 
the 10% mean average percentage error bounds, so that the 
errors of the forecast curve in Fig. [3] are within 10% of the 
actual load. For simplicity of comparison, we do not explicitly 
model DG in this simulation. But we want to emphasize that 
the prediction error will grow larger as the DG penetration 
level increases. Finally, the number of households connected 
to each bus will depend on the test feeder specification for 
each simulation scenario. 
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Base Load Profile 



30% Penetration, IEEE 37 Bus 




11 13 15 17 19 21 23 1 3 5 7 9 11 13 15 
Time (h) 



Fig. 3. Base load profile used in the simulation with IEEE 37 bus system. 



TABLE I 




Vehicle Facts 




Parameter 


Value 


Battery Capacity 


16 kWh 


Energy Usage per 100 miles 


34 kWh 


Charging Rate (120 V, 16 A) 


1.92 kW 


Average Daily Commute Distance 


25 miles 


Daily Consumption 


8.75 kWh 


Charging Efficiency 


0.90 



2) PEV Specifications: The PEV specifications and driving 
patterns are summarized in Table U which are obtained from 
typical PEV specifications l29l and the national transportation 
survey l30ll . The daily consumption of 8.75kWh in Table U is 
obtained from the 25 miles average daily commute distance 
and the PEV consumption rate of 34 kWh/100 miles. The 
average number of registered vehicle is assumed to be 1.8 per 
household |30|. In the simulation, The PEVs will be randomly 
assigned among all vehicles, according to the PEV penetration 
level. The random PEV driving pattern in this simulation is 
based on the national transportation survey 11301 . A PEV is 
assumed to leave home at 7am, with zero-mean Guassian offset 
with standard deviation of 1 hour. Similarly, a PEV arrives 
home at 5pm, with zero-mean Gaussian offset with standard 
deviation of 2 hours. We focus on residential charging in this 
simulation, and assume that the PEV can only charge at home 
after it connects to the grid on returning home. 

B. IEEE 3 7- Bus System 

The PEV charging is first simulated in the standard IEEE 
37-bus test feeder. In this case, the total number of vehicles 
is 3402. There are three types of smart charing algorithms 
considered in the simulation: 1) A static optimal charging 
algorithm which solves d 1 1 1 . with perfect knowledge of the 
day-ahead values of all randomness; 2) A static suboptimal 
charging algorithm, which solves ( fTTT i using imperfect forecast 
of day-ahead load curve; and 3) On-line charging proposed by 
Algorithm Q] in this paper. 

The charging algorithms are simulated at PEV penetration 
levels of 30% and 50%. For the 30% penetration case, /3 = 






Static Optimal 


□ 


Static Suboptimal 


A 


On-line Decentralized 


V 


Base Load 




Fig. 4. The total system loads with 30% PEV penetration in the IEEE 37-bus 
system. 



50% Penetration, IEEE 37 Bus 



4.5 



I 3.5 



Static Optimal 

□ Static Suboptimal 

A On-line Decentralized 

V Base Load 




Fig. 5. The total system loads with 50% PEV penetration in the IEEE 37-bus 
system. 



0.0205, and C, = 577 for each vehicle, whereas for the 50% 
penetration case, (3 = 0.0161, and Ci = 534 for each vehicle. 
In both cases, convergence can be observed after around 20 
iterations for each time slot computation. The maximum total 
computation time of the on-line algorithm is 0.58 second for a 
24-hour simulation scenario, while 3900 seconds for the static 
optimizations. Note the dramatic computation performance 
improvement for the case of on-line charging. This is due to 
the fact that each charging schedule is computed using only 
current system states, which have much smaller dimension 
than the total state processes. In practice, the time scale of each 
time slot is on the order of minutes. Thus, the computation and 
communication requirement of the on-line charging algorithm 
can be easily satisfied. The results of total loads are shown in 
Fig. H]and Fig. [5] respectively. We have the following remarks. 

1) Valley Filling: One can easily verify effectiveness of 
the minimum load variance formulation (O by observing that, 
in both cases, the static optimal solution achieves a perfectly 
flat total load curve at night. Thus, compared to other smart 
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Fig. 6. The total system load profile and typical battery state profile with 
30% PEV penetration in the IEEE 37-bus system with prioritized customers. 
Note that all customers in the 'uniform' case have the same C\ = 577. 



charging formulations, the ones based on electricity price in 
particular, the minimum load variance formulation can avoid 
the additional 'midnight peak', which, in the extreme case, 
may cause similar grid congestion issues as uncoordinated 
charging. 

2) Optimality: The proposed on-line decentralized PEV 
charging achieves almost the same total load profile as the 
static optimal, even though the former does not need to 
know the driving pattern and loads in advance. This further 
verifies the theoretical result in [3] Thus, we can achieve the 
same performance as the static optimal, with much smaller 
computational overhead. 

3) Robustness: The day-ahead prediction based algorithms 
are vulnerable to the forecast errors. This can be clearly 
observed from Fig. |4] and Fig. [5] where the forecast based 
solutions can not achieve a flat profile in the presence of the 
load forecast error. In fact, we allowed these algorithms to 
know the exact driving patterns in advance, which is clearly 
unrealistic. On the other hand, the optimal decentralized charg- 
ing algorithm is not affected by such forecast errors, since it 
is an on-line algorithm, which does not rely on forecasts. 

4) Service Differentiation: The constants {Ci} are used to 
provide service differentiation among the PEVs, where a larger 
Ci implies a higher priority for PEV i. That is, for two PEVs i 
and j with the same SoC, if C; > Cj, PEV i can start charging 
sooner than PEV j, according to cTT~3T > . In order to demonstrate 
this effect, we simulate the 30% penetration case with two 
classes of customers, where 10% customers have high priority, 
with Ci = 877, and 90% customers have low priority, with 
Ci = 577. The total loads and battery states are shown in Fig. 
[6] One can clearly see that the energy queues of both classes 
are 'stable'. Further, the high priority customers can finish 
charging much earlier than the customers with low priorities. 
Note that the values of {Ci] have to be chosen carefully by the 
utility based on system specification and customers' choice, 
which is out of the scope of this paper. 
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Fig. 7. Sensitivity results of the total load with respect to changes in in 
the IEEE 37-bus system. 



5) Effects of (3: The parameter /3 specifies the trade-off 
between charging the PEV batteries, and achieving a small 
charging cost. To illustrate this, we simulate the 30% pene- 
tration case in Fig. □ with /3 = 0.0205,0.0041, and 0.1025, 
respectively. One can clearly observe that, when (3 is small, the 
PEV charging is more greedy, as the weight of the charging 
cost is smaller. On the other hand, large j3 will force the total 
load profile to be more flat, by penalizing the PEV charging 
power. In practice, j3 has to be chosen carefully by the utility, 
using system specification and historical data. 

C. IEEE 123-Bus System 

We next simulate the PEV charging in IEEE 123-bus test 
feeder, where the total number of vehicles is 4832. In this 
case, one cannot implement the static optimizations on the 
machines specified at the beginning of this section, due to 
the large problem size (around 4600 thousand variables). On 
the other hand, the on-line charging can still be implemented, 
as the charging decisions are made only using the current 
system states, which have much lower dimension (around 4.8 
thousand variables). The charging algorithms are simulated 
at PEV penetration levels of 30% and 50%. For the 30% 
case, 13 = 0.0205, d = 815, while for the 50% case, 
(3 = 0.161, = 764. In both simulations, convergence can 
be reached after around 20 iterations. The total computation 
time of the on-line algorithm is 0.75 second for a 24-hour 
simulation scenario. The results of total loads are shown in 
Fig. [8] and Fig. [9] respectively. One can easily observe that, 
in both cases, the on-line algorithm can achieve a flat load 
profile. Thus, the on-line algorithm not only achieves the 
optimal charging cost, but also is dramatically computationally 
efficient, and much easier for implementation in large-scale 
systems, compared to other approaches forecast and static 
optimization based approaches. 

V. Conclusion 

This paper proposed an on-line decentralized PEV charging 
algorithm to minimize the distribution system total load vari- 
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Fig. 8. The total system load under the on-line decentralized algorithm with 
30% PEV penetration in the IEEE 123-bus system. 
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Fig. 9. The total system load under the on-line decentralized algorithm with 
50% PEV penetration in the IEEE 123-bus system. 



ance. The charging decisions are made locally at each PEV 
based entirely on the power system states in each time slot. 
Since driving pattern and load forecasts are not needed in this 
algorithm, the charging performance is robust against various 
uncertainties in the system, as compared to the prediction 
based static optimization approaches in the literature. In the 
analysis, asymptotic results about the cost function and battery 
level stability are shown. The performance of the proposed 
charging algorithm is further compared against static smart 
charging algorithms by simulation results. 



Appendix A 
Proof of TheoremQ] 

Proof: It is obvious that any feasible solution of (fTTT i is 
also feasible for (O-®. Now, it is sufficient to show that the 
two objective functions are different by at most a constant that 
is independent of choice of charging schedules {Pj(n)}. Note 



that the constraint © implies that 



(19) 



n=l 



This, together with the definition of /i in (TlOb imply that 



1 T N 



(20) 



Therefore, we have 
2// 



i=i 



N 



/ = /-^£(Sne t (n)+£p(n))+ M 2 (21) 



T 

5 



i=l 



(22) 



where, by ( flOl ), fi is independent of the choice of schedules. 
Thus, the theorem holds. ■ 

Appendix B 
Proof of Theorem[2] 

Proof: We can formulate the optimization in ( fl3] l as the 
following equivalent problem: 

N 

max Y(Ui(n- I) + a) Vi P t At - f3y 2 (23) 



N 



subject to y = ^ P + 5 net (n) 



(24) 

(25) 
(26) 
(27) 



Q < Pi < Pi max , 1 < i < N 
PiT/.At < Ui(n - 1), 1 < i < N 
Pi = if Xi(n) = 0, 1 < i < N 

Thus, one can easily write the dual problem as 

N 

min max ^({7,(n - 1) + C^^PjAt - (3y 2 
Umi {Pi,y} r— f 

JV 

+C/ref(2/-^P-5 net (n)) 
i=l 

subject to l|25]l. (|26|) and ||27j (28) 

It can be easily shown that y* = U Te f/2(3. Thus, we can 
simplify (l28l as 



N 

min max > (UAn — 1) + C,; ^— )riiPiAt 



subject to (|25|), |26| and (|27| 



Kef 

(£/ref) 2 

4/3 



- Kef ^net (™) 



(29) 



One can immediately see that the decentralized charging in 
( TT3l l corresponds to the inner optimization problem with charg- 
ing power constraint d25t . and the update on U le f in Step 3 of 
Algorithm Q] corresponds to standard binary search algorithm 
over the scalar U le f to solve the convex dual problem. ■ 
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Appendix C 
Proof of Theorem[3] 

In this section we prove Theorem [3] Define a 'Lyapunov' 
type function F(n) as follows: 

N n N 

F ( n ) = 2 E(w + C ^ + /'Eowao + E p ^ k )f 

i=l k=l i = l 

The key in proving both cost optimality as well as stability 
results involves analyzing the drift behavior of F(n), which 
we will illustrate in the following. Note that similar techniques 
have also been applied in [16] in the context of wireless 
networks. 

A. Cost Optimality 

The proof of ( [T7| > requires several technical lemmas. The 
following lemma provided a bound on the single slot drift of 

F(n). 

Lemma 1: The one-slot drift of F(n) can be bounded as 
follows: 

AxF(O) 
= F(l)-F(Q) 

N 

< 53(tfi(0) + Ci){Ai{l) - %Pi(l)At) + B 3 
i=i 

N 

+/3(Sne,(l)+^P(l)) 2 (30) 



Proof: For each PEV i, direct calculation shows that 

1(11,(1) + ar 
= ^(Ui(o) + a ( (i) - 77^(1) a* + a) 2 

= ~(Ui(fl) + C t ) 2 + (Ui(0) + C l )(A l (l) - rnWAt) 



+-(A l (l)-r 1% P l (l)At) 2 (31) 
Thus, the following bound holds: 

+ Ci ) 2 - 1(^(0) + cy 2 

= (Ui(p) + Ct)(AiQ-)-ViPiO)*t) 

+±(A i (l)-r H P i (l)At) 2 

< (Ui(0) + -7,^(1) At) + ^Kf (32) 

where K t = max(A™ ax , i]iP™ & *At). The last inequality holds 
because we assume that a PEV can either charge (in garage) or 
discharge (on the road) in each time slot, but not both. Thus, 
the lemma follows after summing up the above inequality 
over N customers and adding the charging cost drift term 

We next generalize the above bound to T time slots. 
Lemma 2: Under Algorithm[T] the first T-slot drift of F(ri) 
can be bounded as 

T(T+1) 



Proof: Carrying out the procedure in Lemma Q] over T 
time slots, we obtain the following 



T N 



A T F(0) < ^^(Uiin-y + CiXAiinj-rjiPiWAt) 

n—l i—1 

T N 

+B 3 T + /?^(S net (n)+^P(n)) 2 (34) 

n—l i—1 

Now, we want to compare the bound on the RHS of d34l ) under 
the myopic charging schedule {Pi(n)}, which are computed 
by Algorithm Q] against a static optimal solution of (fTTT i 
with perfect knowledge of all stochastic dynamics, which we 
denote as {i^ stat (n)}. Because Algorithm Q] always greedily 
maximizes ( fTBT l, the RHS of the above can be further bounded 
as follows: 

A T F(0) 

T N 

< E E^ - X ) + C *)(M«>) - ViPr(n)At) 

n—l i—1 

T N 

(B 1 + B 3 )T + ^(S net (n) + J2 A smt («)) 2 (35) 

n=l i=l 

The extra term B\ = YjiLiiQ + Ki)K{ is introduced to 
bound the case where the energy queues are small, so that the 
charging process stops in the middle of a time slot. In such a 
case, we have 

{Ui(n ~l) + d)THP?*(n)At < (K t + C t )Ki (36) 

We next bound the drift in the energy queue lengths. Note that 
the energy queue length for each PEV i can be written as the 
following: 

n 

Ui(n) = 11,(0) + (Mn) ~ ViPi{n)At) (37) 

k=l 

Thus, Ui(n) can be bounded as follows: 

U^O) - K in < U^n) < U^O) + K in (38) 

which yields the following bound: 

Ui(n-l)(Ai(n)-rhPrin)M) < 

Ui(0)(Ai(n) - ViPnn)At) + K 2 (n - 1) (39) 

Plugging the above inequality into the d35l l yields 

A T F(0) 

N T 

< E(^(°) + ^(E (Mn) ~ ^Pr(n)At)) + B\T 



i=i 



N 



-B 3 T + B 3 £(n - 1) + ^(SU(«) + E P ">^ 



A T F(0) < B X T + B 3 



+ 0f*(l)T (33) 



® Bl T + B 3 T{T n +1) + 0f*(l)T 



(40) 
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where step (a) is because {P* tat (n)} solves the optimization 
(fTTT i. and satisfies 



(Mn) ~ wifHAt) = (41) 

T JV 



follows: 



Thus, the lemma follows. ■ 
We next extend the above analysis from one day to D days: 

Lemma 3: The drift of F(n) over the first DT time slots 
satisfies 

A DT F(0) <B 1 DT + B 3 ^^-^-+l3Tj2f*(d) 



Proof: This can be simply obtained by summing the 
bound in Lemma |2] over D days. ■ 
We are now ready to prove the first part of Theorem [3] 

Proof: (Cost Optimality of Algorithm^ According to the 
bound in Lemma [3] the average cost of the first D days under 
Algorithm Q] can be bounded as 



DT N 



DT 



n=l 



< F(DT) 

~ (3DT y ' 

_ F(0) + A DT F(0) 

(3DT 



D 



< 



*W>.* +a r+i + »Er M (43) 



DT p 



2/3 D 



Thus, after taking D — > oo, the cost optimality claim in ( fTTI i 
follows. ■ 



B. Stability 

The proof of the second part of Theorem[3]requires another 
set of lemmas bounding the drift of F(n). We start with the 
first one the T-slot drift of the function F{n). 

Lemma 4: Under AlgorithmQ] the first T-slot drift of F(n) 
can be bounded as 

T N 

A T F{0) < -e^2^2(U i (n)+C i )+B 1 T 

n—l i—1 

T(T + 1) 

+(B 2 + B 3 ) 1 ^ + /3/ max T (44) 



Proof: Denote {-ff as (n)} as a sequence of strictly fea- 
sible charging schedules, which satisfy ([Tol l. Because the 
charging schedules {Pj(n)} by Algorithm Q] always greedily 
maximize Sl5[ , we can obtain a similar bound as d35l l as 



A T F(0) 

T N 

< E E^( n - !) + G 0(A(n) - rf a »At) 



n=l (=1 



AT 



+BiT + B 3 T + /3]T(S net (n) + E ^M) 2 



T iV 



< E E(^( n - x ) + ^)(^(") - rf a »At) 

n—l i—1 

+B X T + B 3 T + ^/ max T (45) 

From d38l l. the following bound holds: 

Ui(n - l)(Mn) - Vl P^(n)At) < 
Ui(Q)(Ai(n) ~ ^ ea »At) + K?(n - 1) (46) 

Plugging the above into (l45l l and applying (fTol l yields 

A T F(0) 

N T 
i—1 n—l 

T{T + 1) 



-B{T + B 3 



/3/ max r 



< -sTj2(Ui(0) + d) + B 3 T[T + l) 



+B{T + pf^T (47) 
Finally, from (l38l the above can be further bounded as 

A T F(0) 

N T 
i—1 n—l 

T(T+1) 
+BiT + B 3 2 + /3/ max T 

i=l n—l 

T(T + 1) 

+BiT + B 3 2 + /3/ max T (48) 

from which the lemma follows. ■ 
Similarly, the drift of F(n) over D days can be bounded 

by the following lemma. 

Lemma 5: The drift of F(n) over the first D days can be 

bounded as follows: 

DT N 

A DT F{0) < -e^^ U ^ n ) + C ^ +BlDT 

n—l i—1 

+ (B 2 + B 3 ) DT{T + l) + pDf m ^T (49) 



Proof: This can be easily verified by applying Lemma [4] 
D times and summing the bounds. ■ 
We are now ready to prove the second part of Theorem [3] 
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Proof: (Stability of Algorithm^ From d49b one can easily 
see that 



< 



< 



DT N 
n— 1 i—1 

F(0)~F(DT) B 1 (B 2 + B 3 )(T+l) t3f 

eDT e 2e 



£ 



F(0) , B x , (B 2 + B 3 )(T + 1) , 



(50) 



from which (fT8l l follows from above by taking D — > oo. ■ 
References 

[1] A. Boulanger, A. Chu, S. Maxx, and D. Waltz, "Vehicle electrification: 
Status and issues," Proc. IEEE, vol. 99, no. 6, pp. 1116-1138, Jun. 2011. 

[2] G. Heydt, "The impact of electric vehicle deployment on load manage- 
ment strategies," IEEE Trans. Power App. Syst., vol. PAS-102, no. 5, 
pp. 1253-1259, May 1983. 

[3] S. Meliopoulos, J. Meisel, G. Cokkinides, and T. Overbye, "Power 
system level impacts of plug-in hybrid vehicles," Power Systems En- 
gineering Research Center (PSERC), Tech. Rep., Oct. 2009. 

[4] K. Clement-Nyns, E. Haesen, and J. Driesen, "The impact of charging 
plug-in hybrid electric vehicles on a residential distribution grid," IEEE 
Trans. Power Syst., vol. 25, no. 1, pp. 371-380, Feb. 2010. 

[5] L. Fernandez, T. Roman, R. Cossent, C. Domingo, and P. Frias, 
"Assessment of the impact of plug-in electric vehicles on distribution 
networks," IEEE Trans. Power Syst, vol. 26, no. 1, pp. 206-213, Feb. 
2011. 

[6] J. Lopes, F. Soares, and P. Almeida, "Integration of electric vehicles in 
the electric power system," Proc. IEEE, vol. 99, no. 1, pp. 168-183, 
Jan. 2011. 

[7] K. Dyke, N. Schofield, and M. Barnes, "The impact of transport elec- 
trification on electrical networks," IEEE Trans. Ind. Electron., vol. 57, 
no. 12, pp. 3917-3926, Dec. 2010. 

[8] K. Schneider, C. Gerkensmeyer, M. Kintner-Meyer, and R. Fletcher, 
"Impact assessment of plug-in hybrid vehicles on pacific northwest 
distribution systems," in Proc. Power Energy Soc. Gen. Meet, Jul. 2008, 
pp. 1-6. 

[9] D. Wu, D. Aliprantis, and K. Gkritza, "Electric energy and power 
consumption by light-duty plug-in electric vehicles," IEEE Trans. Power 
Syst, vol. 26, no. 2, pp. 738-746, May 2011. 

[10] L. Gan, UT, and S. Low, "Optimal decentralized protocols for electric 
vehicle charging," in Proc. IEEE Conf. on Dec. and Contr, 2011. 

[11] Z. Ma, D. Callaway, and I. Hiskens, "Decentralized charging control for 
large populations of plug-in electric vehicles: Application of the nash 
certainty equivalence principle," in Proc. IEEE Conf. on Dec. and Contr. , 
Sep. 2010, pp. 191-195. 

[12] E. Sortomme, M. Hindi, S. MacPherson, and S. Venkata, "Coordinated 
charging of plug-in hybrid electric vehicles to minimize distribution 
system losses," IEEE Trans. Smart Grid, vol. 2, no. 1, pp. 198-205, 
Mar. 2011. 

[13] N. Rotering and M. Ilic, "Optimal charge control of plug-in hybrid 
electric vehicles in deregulated electricity markets," IEEE Trans. Power 
Syst, vol. 26, no. 3, pp. 1021-1029, Aug. 2011. 

[14] D. Bertsekas, Dynamic programming and optimal control. Athena 
Scientific, 1995. 

[15] S. Meyn, Control techniques for complex networks. Cambridge 
University Press, 2008. 

[16] M. Neely, Stochastic network optimization with application to commu- 
nication and queueing systems. Morgan Claypool, 2010. 

[17] N. McKeown, A. Mekkittikul, V. Ananfharam, and J. Walrand, "Achiev- 
ing lOOthroughput in an input-queued switch," IEEE Trans. Commun., 
vol. 47, no. 8, pp. 1260-1267, Aug. 1999. 

[18] Q. Li and R. Negi, "Scheduling in wireless networks under uncertainties: 
A greedy primal-dual approach," in Proc. IEEE Int. Conf. Commun. , Jun. 
2011, pp. 1-5. 

[19] , "Back-pressure routing and optimal scheduling in wireless broad- 
cast networks," in Proc. IEEE Global Commun. Conf, Dec. 2009, pp. 
1-6. 

[20] A. Saber and G. Venayagamoorthy, "Plug-in vehicles and renewable 
energy sources for cost and emission reductions," IEEE Trans. Ind. 
Electron., vol. 58, no. 4, pp. 1229-1238, Apr. 2011. 



[21] O. Sundstroem and C. Binding, "Optimization methods to plan the 
charging of electric vehicle fleets," in Proc. Int. Conf. on Control, 
Commun. and Power Eng., 2010. 

[22] S. Deilami, A. Masoum, P. Moses, and M. Masoum, "Real-time coor- 
dination of plug-in electric vehicle charging in smart grids to minimize 
power losses and improve voltage profile," IEEE Trans. Smart Grid, 
vol. 2, no. 3, pp. 456-467, Sep. 2011. 

[23] K. Turitsyn, N. Sinitsyn, S. Backhaus, and M. Cherfkov, "Robust 
broadcast-communication control of electric vehicle charging," in Proc. 
IEEE SmartGridComm, Oct. 2010, pp. 203-207. 

[24] R. Lasseter, "Smart distribution: Coupled microgrids," Proc. IEEE, 
vol. 99, no. 6, pp. 1074-1082, Jun. 2011. 

[25] S. Han, S. Han, and K. Sezaki, "Development of an optimal vehicle- 
to-grid aggregator for frequency regulation," IEEE Trans. Smart Grid, 
vol. 1, no. 1, pp. 65-72, Jun. 2010. 

[26] W. Kersting, "Radial distribution test feeders," in Proc. IEEE Power 
Eng. Soc. Win. Meet, vol. 2, 2001, pp. 908-912. 

[27] S. C. Edision. (2011) Regulatory information - see load profiles. [On- 
line]. Available: http://www.sce.com/AboutSCE/Regulatory/loadprofiles/ 

[28] U. E. I. Administration. (2009) Residential average monthly 
bill by census division and state. [Online]. Available: 
http://www.eia.gov/cneaf/electricity/esr/table5.html 

[29] U. E. P. Agency. (2011) 2011 electric vehicles, fuel economy. [Online]. 
Available: http://www.fueleconomy.gov/feg/evsbs.shtml 

[30] P. Hu and T. Reuscher. (2004, Dec.) Summary of travel trends. U.S. 
Department of Transportation and Federal Highway Administration. 
[Online]. Available: |http://rmts.oml.gov/20 01/pub/STT.pdf] 



Qiao Li (S'07) received the B.Engg. degree from the Department of Elec- 
tronics Information Engineering, Tsinghua University, Beijing China, in 2006. 
He received the M.S. degree from the Department of Electrical and Computer 
Engineering, Carnegie Mellon University, Pittsburgh, PA USA, in 2008. He 
is currently a Ph.D. Candidate in the Department of Electrical and Computer 
Engineering, Carnegie Mellon University. His research interests include dis- 
tributed algorithms, smart grids technologies, and wireless networking. 



Tao Cui (S'10) was born in Shaanxi Province, China, in 1985. He received 
the B.Sc. and M.Sc. degrees from Tsinghua University, Beijing, China, and is 
currently pursuing the Ph.D. degree at Carnegie Mellon University, Pittsburgh, 
PA. His main research interests include power system computation, protection, 
analysis, and control. 



Rohit Negi received the B.Tech. degree in electrical engineering from the 
Indian Institute of Technology, Bombay, in 1995. He received the M.S. and 
Ph.D. degrees from Stanford University, CA, in 1996 and 2000, respectively, 
both in electrical engineering. Since 2000, he has been with the Electrical and 
Computer Engineering Department, Carnegie Mellon University, Pittsburgh, 
PA, where he is a Professor. His research interests include signal processing, 
coding for communications systems, information theory, networking, cross- 
layer optimization, and sensor networks. Dr. Negi received the President of 
India Gold Medal in 1995. 



Franz Franchetti received the Dipl.-Ing. degree and the PhD degree in 
technical mathematics from the Vienna University of Technology in 2000 
and 2003, respectively. Dr. Franchetti has been with the Vienna University 
of Technology since 1997. He is currently an Assistant Research Professor 
with the Dept. of Electrical and Computer Engineering at Carnegie Mellon 
University. His research interests concentrate on the development of high 
performance DSP algorithms. 



12 



Marija D. Hie (M'80-SM'86-F'99) is currently a Professor at Carnegie Mel- 
lon University, Pittsburgh, PA, with a joint appointment in the Electrical and 
Computer Engineering and Engineering and Public Policy Departments. She is 
also the Honorary Chaired Professor for Control of Future Electricity Network 
Operations at Delft University of Technology in Delft, The Netherlands. She 
was an assistant professor at Cornell University, Ithaca, NY, and tenured 
Associate Professor at the University of Illinois at Urbana-Champaign. She 
was then a Senior Research Scientist in Department of Electrical Engineering 
and Computer Science, Massachusetts Institute of Technology, Cambridge, 
from 1987 to 2002. She has 30 years of experience in teaching and research 
in the area of electrical power system modeling and control. Her main interest 
is in the systems aspects of operations, planning, and economics of the electric 
power industry. She has co-authored several books in her field of interest. Prof. 
Ilic is an IEEE Distinguished Lecturer. 



